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Abstract. - The nonadiabatic regime of the electron-phonon interaction leads to behaviors 
of some physical measurable quantities qualitatively different from those expected from the 
Migdal-Eliashberg theory. Here we identify in the Pauli paramagnetic susceptibility x of 
such quantities and show that the nonadiabatic corrections reduce x with respect to its adiabatic 
limit. We show also that the nonadiabatic regime induces an isotope dependence of x, which in 
principle could be measured. 



When the Fermi energy Ep is anomalously small, as in high- Tc cuprates ]l|] and in the 
fuUerene compounds the Migdal-Eliashberg (ME) approach ^ may result inadeguate in 
describing the interplay between charge carriers and phonons. For example, the alkali-doped 
fuUerenes (A3C60) have Fermi energies of order 0.25 eV Q and intramolecular phonon modes 
with frequencies wo in the range between 20 meV and 0.2 eV 1^. In this case, the adiabatic 
parameter loq / Ep lies somewhere between 0.1 and 0.9, depending on which phonon modes most 
couple to the electrons. The main consequence is that the electron-phonon vertex corrections 
may no longer be negligible, as assumed in the ME framework, and a generalization of the 
theory is required to include the nonadiabatic contributions 

This generalization In terms of the electron-phonon coupling A and the adiabatic parameter 
cuo/Ep, the ME regime applies for A <, 1 and loq/Ep <C 1. Therefore, a generalization beyond 
the ME framework is required when A> 1 and/or ljq/Ep is no longer negligible. However, 
when A is larger than some critical value Ac (which is of order one or larger), the system 
evolves toward a polaronic regime characterized by strong electron-lattice correlations. This 
holds true even in the adiabatic case in which the charge carriers aquire large effective masses. 
On the other hand, a region in the X-Uo/Ep plane different from the one leading to polaronic 
states is defined by A ;< 1 and LOg/Ep finite. Within this region, where the charge carriers 
are weakly interacting nonadiabatically with phonons, the nature of quasiparticles is different 
from both the ME and the polaronic ones. In such a nonadiabatic regime we shall speak of 
nonadiabatic Fermi liquids (or nonadiabatic fermions), to stress the difference from the ME and 
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the polaronic limits. In practice, such a regime can be described by a perturbative approach 
where Xloq/Ef plays the role of the small parameter of the theory j6[ |^. Various comparisons 
with exact results (for the one electron case) ^ and quantum Monte Carlo calculations [|| 
point toward the reliability of such a perturbative description. 

At the zeroth order in Xloq/Ef, the nonadiabatic theory coincides with the ME limit while 
for finite values of Xluq/Ep the nonadiabatic fermions display anomalous behaviors. In this 
situation, several properties are modified and a very important question regards the possibility 
to observe some fingerprints of such a nonadiabatic regime. Furthermore, in order to be 
considered as possible evidences, such fingerprints should be searched among those physical 
quantities for which some well established property in the ME regime results qualitatively 
modified in the nonadiabatic one. In order to clarify this statement, let us consider for example 
the electron-phonon renormalized charge carrier mass m*. In the ME regime m* = (1 + A)m 
, where m is the bare mass and A is the electron-phonon coupling. Since A is independent 
of the ion-mass no isotope effect is expected for m*. However, when the nonadiabatic 
contributions are no longer negligible, m* aquires an ion-mass dependence which leads to 
a non-zero isotope coefficient am* (l^- The effective mass m* represents therefore a clear 
example of a quantity for which a well established property in the ME regime (a™* = 0) is 
drastically modified in the nonadiabatic one (am* 7^ 0). So far, strong evidences for isotope- 
dependent m* have been reported for YBa2Cu306+2: and La2_xSr2,Cu04 and theoretical 
calculations have shown that already the inclusion of the first nonadiabatic vertex correction 
to the ME limit provides values of with sign and order of magnitude in agreement with 
those estimated by the experiments . 

Another property typical in the ME regime which is instead strongly altered by the nona- 
diabatic contributions is the non-magnetic impurity dependence of the critical temperature Tc 
of an homogeneous s-wave superconductor. For a conventional superconductor, weak disorder 
does not influence the critical temperature as stated by Anderson's theorem |Q. On the 
contrary, since the electron-phonon vertex corrections are very sensitive to the amount of 
disorder, the critical temperature of a s-wave nonadiabatic superconductor can be strongly 
lowered by the impurities fl^ j . Such a peculiar behavior is also accompained by an anomalous 
impurity dependence of the isotope coefficient of T^. So far, reduction of TJ, driven by disorder 
for s-wave superconductors has been reported for K^Cqq @ and Nd2_a;Ce3Cu04_5 

In this paper we consider another measurable quantity which could be considered as a test 
for the breakdown of Migdal's theorem: the Pauli paramagnetic susceptibility x- Here, the 
characteristic feature in the ME regime {<jJq/Ef <^ 1) is that the electron-phonon interaction 
does not renormalize the Pauli susceptibility so that x is independent of A and cjq [|io|. In 
the ME regime therefore x = Xp — where is the Bohr magneton and A^(0) 

is the electron density of states at the Fermi level. In principle, therefore, a measure of x 
via for example electron paramagnetic resonace (EPR) is unaffected by the electron-phonon 
interaction and provides an estimate of the electronic density of states iV(0) which however is 
renormalized by many-electrons effects (Stoner enhancement) (EJ). 

The interesting aspect of x is that, as we show below, when Lo^jEp is no longer negligible x 
aquires a phonon renormalization and becomes dependent on both A and loq. This result can 
be of importance for two reasons. First, it leads to re-consider the estimates of the electron 
density of states obtained by EPR measurements, since these estimates have been based on 
the phonon-independent ME form of x- Second, and more importantly, the nonadiabatic 
renormalization of x induces a non-zero isotope effect which, in principle, could be measured. 

To evaluate the Pauli susceptibility we make use of the static limit of the Kubo formula: pSl 



(^) In the present discussion we shall consider the many-electrons effects as being already contained 
in Af(0). 
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X(r)= lim ul / dT{TrSM,r)S,{-q,0)), 
Jo 

where P is the inverse temperature T and 



(1) 



SM= E -U.^k.. (2) 

k,CT=±i 

where ("^ktr) creation (annhiUcation) operator for electron with momentum k and 

spin direction cr = ±1. 

In what foUows, we shaU focus on the evaluation of Eq.(|l]) for a system of electrons 
interacting with phonons through the coupling (?(q). In terms of electron and phonon Green's 
functions, Eq.(|l|) reduces to the following general expression: 

X(T) = - hm mI^E E k + ^^^^^ + "^)' (3) 

m k 

where a;,„ = (2m + 1)ttT and 

G(m,k) = [ic^„-e(k)-S(m,k)]-\ (4) 

is the Green's function for an electron with dispersion e(k) and electron-phonon self-energy 
S(m. k). In Eq.(^), r(k + q, k; m) is the irreducible electron-phonon vertex function which is 
given by all diagrams which cannot be separated into two different parts by cutting a single 
electron or phonon propagator line. The reducible part of the vertex function gives in fact 
zero contribution when the summation over the spin indeces is performed in eqs.(0j^) ]T^ . 

In this paper we compute eq. by employing a self-consistent calculation which amount to 
evaluate I](m, k) in the non-crossing approximation. For dispersionless phonons with frequency 
Wo, we consider therefore the electron-phonon self-energy as given by: 



S(r., k) = T E 3(k - k')^7 ^G{m, k'). (5) 

mk' 

In the above equation we have implicitly assumed that the phonons are already renormalized 
and that loq is a dressed phonon frequency. In a conserving approach, the vertex function 
resulting from the non-crossing approximation for E(m, k) is given by all the ladder contribu- 
tions. Therefore the vertex function satisfies the following ladder equation: 



r(k-^q,k;n-^TO,n) = 1+ tV .g(k-k')^7 ^ ^Gfm', k')G(m' + m, k' + q) 

, ■ (t^n - t^m'j + ^0 

m'K' 

xr(k'-Hq,k';m' + TO,TO'). (6) 

Actually, from eq.(^, to evaluate x we only need to retain the static limit of eq.(|^) which is 
given by setting first = and after q = 0. As already shown in Refs.[0, if we exchange the 
order of the two limits, the resulting dynamical limit of the vertex will be in general different 
from the static one. Therefore setting LUm = and q = in both hand sides of eq.(P) may give 
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a non well defined result because in that point the vertex in non-analytic. However, as we shall 
show below, the computing procedure we employ in handling the vertex function automatically 
provides the correct static limit by simply setting LOm = 0, q = in eq.(^, regardless of the 
order of the two limits. Therefore, by setting limq^or(k + q, k; n, n) ~ rs(k, 7i), the static 
limit of eq.(|^) reduces to: 



r,(k,n)=.l + TV .g(k-k')27 ^^t^-^— ^G(m', k')2r,(k', m'). (7) 

m'k' 

Without loss of generality, the solution of the set of equations (^) , (||) and (0) can be found 
by using a structureless electron-phonon interaction .g(q) = . The resulting self-energy is 
then momentum independent and for a system with a half-filled electron band of constant 
DOS over the entire bandwidth 2Ep, the self-energy can be written as = iw„ — iW^, 

where 



Wn = ujn + AvrT y „ , 2 - arctan ( ^ ) , (8) 

is the renormalized electron frequency and A = g^N{0) is the electron-phonon coupling. Within 
the same approximation scheme, rs(k, n) becomes momentum-independent and the resulting 
vertex function Ts{n) satisfies the following equation: 



2 f Ef 



We can verify that the above equation gives indeed the static limit of the vertex by neglecting 
the renormalization of the frequency, W^' '^m', and by performing the zero temperature 
limit. In this way, to the first order in A and at zero external frequency, Eq.(0) becomes: 



r.(0) = 1 - A / ^^^^^ = 1 - A— (10) 

rs(0) coincides therefore with the static limit already calculated in the perturbation theory. 
I 

We are now in the position to evaluate the Pauli susceptibility. Since both the self-energy 
and the vertex function are independent of the momentum, equation can be analitically 
integrated over the energy and the final expression for x(r) reduces to: 



xm-^'^El^r.M, (11) 

where xp ~ Mb-^(0)i and Wm and Tsim) are the solution of equations (||) and (^), respectively. 

We solve the set of equations (H), (|^) and ( pT| ) for a temperature T/loq = 0.02 and differtent 
values of A and loq/Ef- The frequency summations appearing both in the self-energy (|^) 
and in the vertex function (^ is truncated at the frequency cut-off = {2N + 1)'kT with 
N = 400 corresponding to ujc — 50ujq. The solutions of (|^) and (^) are then calculated by 
iteration and the results are plugged into eq. (pT|) . The high-frequency part {um > cUc ^ ^o) 
of the summation in eq.(0) is calculated by setting Wm = i^m and Ts{m) = 1, since in this 
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high-frequency region the contribution from the electron-phonon couphng is neghgible. The 
procedure outhned above permits to estimate the zero temperature susceptibihty x s-lso for 
the smallest value of ujq/Ef we used in the calculations {ojq/Ef = 0.01). 

In Fig. 1 we show the zero temperature calculated Pauli susceptibility as a function of the 
adiabatic parameter ujq/Ep and for different values of the electron-phonon coupling constant 
A. When ujq/Ep — > 0, x approaches its free-electron value xpi irrespectively of the value of 
A and we recover therefore the result of the ME theory. Instead, when loq/Ep is larger than 
zero, X becomes dependent of A and results to be always lowered with respect to xp- I^i Fig. 2 
x/Xp is plotted as a function of the electron-phonon coupling A for different values oiujo/Ep. 
For small values of ujq/Ep, x/xp decreases almost linearly with A. The main result of our 
calculations is therefore that x(0)/xp < 1 as soon as ujo/Ep > 0. Preliminary calculations 
including higher orders vertex corrections confirm this feature. 

The reduction of the Pauli susceptibility induced by the electron-phonon interaction when 
loq / Ep is finite requires to re-consider the estimates of the electron density of states based on 
EPR measurements In these estimates, in fact, the measured x is fitted with the ME 

expression of the susceptibility 



XC.A'W-^i^, (12) 

where in the last equality we have explicitly separated -/V(O) into the free-electron form A^o(O) 
and the Stoner enhancement 1/(1 — /) given by the many- electrons effects. Theoretical esti- 
mations of 1/(1 — /) permit therefore to obtain A'o(O) from the experimental x However, 
this procedure may sistematically underestimate N^ifi) if ljq / Ep is no longer negligible like in 
the fullerene compounds. In fact, in view of the previous results, A^o(O) of eq.([12|) should be 
replaced by N^{{)) ~ A^o(0)/(A,tJo/-E'F), where the function / takes into account the phonon 
renormalization effects and is less than the unity. From the calculations showed in Figs. 1 
and 2, / can be as small as ^ 0.8 — 0.7, leading to an underestimation of the bare density of 
states 7Vo(0) of - 20 - 30%. 

Another remarkable feature of the nonadiabatic phonon renormalization is the lattice in- 
duced isotope effect on x- From Fig.l in fact it is obvious that a change in frequency loq 
induces a lowering of x- Such a change of loq can be induced by isotope substitution leading 
therefore to a non-zero value of the isotope coefficient: 



^ dlogx ^ 1 dlogx , . 

"'^ dlogM 2d\og{u;o/Ep)' ^ ' 

where M is the ion mass and, in the last equality, we have used ujq cx (Af)~^/^ (note that 
in the nonadiabatic regime x depends also on A, however A is independent of M). In Fig. 3 
we show the numerical evaluation of eq.|l^ as a function of ojq/Ep and for different values of 
A. As expected, the resulting isotope coefficient vanishes at the adiabatic limit. However, 
for nonzero values of Lo^jEp^ it becomes negative and for ordinary values of A can be of 
order —0.05. This is a rather small value, nevertheless it provides a clear indication of 
nonadiabaticity. It would be extremely interesting to investigate experimentally the presence 
or the absence of an isotope effect on x in the fullerene compounds. The outcome of such 
kind of experiment could provide us with an estimate of uq/ Ep and therefore of the degree of 
nonadiabaticity in such narrow band materials. 
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Fig. 1. uoq/Ep dependence of the Pauli susceptibility % for different values of the electron- 
phonon coupling constant A. 
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Fig. 2. - Pauli susceptibility x ^is a function of the electron-phonon coupling constant A for 
different values of the adiabatic parameter loq/Ef- 




Fig. 3. - Isotope cocfRcient of the Pauli susceptibility for different values of the electron- 
phonon coupling A. 



